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Abstract 



We present a new certified and complete algorithm to compute arrangements of real planar algebraic 
curves. Our algorithm provides a geometric-topological analysis of the decomposition of the plane 
induced by a finite number of algebraic curves in terms of a cylindrical algebraic decomposition of the 
plane. Compared to previous approaches, we improve in two main aspects: Firstly, we significantly 
reduce the amount of exact operations, that is, our algorithms only uses resultant and gcd as purely 
symbolic operations. Secondly, we introduce a new hybrid method in the lifting step of our algorithm 
which combines the usage of a certified numerical complex root solver and information derived from the 
resultant computation. Additionally, we never consider any coordinate transformation and the output 
is also given with respect to the initial coordinate system. 

We implemented our algorithm as a prototypical package of the C++-library Cgal. Our implemen- 
tation exploits graphics hardware to expedite the resultant and gcd computation. We also compared 
our implementation with the current reference implementation, that is, Cgal's curve analysis and 
arrangement for algebraic curves. For various series of challenging instances, our experiments show that 
the new implementation outperforms the existing one. 

1 Introduction 

Computing the topology of a planar algebraic curve 



can be considered as one of the fundamental problems in real algebraic geometry with numerous applications 
in computational geometry, computer graphics and computer aided geometric design. Typically, the 
topology of C is given in terms of a planar graph Qq isotopic to C. For a geometric-topological analysis, 
we further require the vertices of Qc to be located on C. In this paper, we study the general problem 
of computing an arrangement of a given set of algebraic curves, that is, the decomposition of the plane 
into cells of dimensions 0, 1 and 2 induced by the given curves. The proposed algorithm is certified and 
complete, and the overall arrangement computation is exclusively carried out in the initial coordinate 
system. Efficiency of our approach is shown by implementing our algorithm and comparing it to the 
current reference implementation. 

There exist a number of certified and complete approaches to determine the topology of an algebraic 
curve; we refer the reader to \12\ \18 \ [231 , 125^ I29j for recent work and further references. At present, only 
the method from |18j has been extended to arrangement computations of arbitrary algebraic curves |17j . 
Common to all existing approaches is that, in a first step, they use eliminations techniques (e.g., resultants) 
to project the x-critical points (i.e., points p £ C with fy{p) = 0) of the curve into one dimension. In 
a second step, the fiber at each of these projected points is computed. In general, this lifting step has 
turned out to be the most time-consuming part because it amounts to determining the real roots of a 
non-square univariate polynomial f{a,y) S M[?/] with algebraic coefficients. The high computational cost 
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for computing the roots of /(a, y) is mainly due to a more comprehensive algebraic machinery such as 
subresultants (in [13 HEl ES] ) , Grobner basis or rational univariate representation (in [12]) in order to 
obtain additional information on the number of distinct real (or complex) roots of f{a, y) or the multiplicity 
of the multiple roots of f{a,y). In addition, all except the method from |12j consider a shearing of the 
curve which guarantees that the sheared curve has no two x-critical points sharing the same x-coordinate 
which in turn simplifies the lifting step but for the price of giving up sparseness of the initial input; an 
approach, which usually results in larger bitsizes of the coefficients and considerably increased running 
times. In a final connection step, arcs adjacent to the same x-critical points are identified. 

The high-level description of the algorithm presented in this paper is almost identical to that of the 
existing methods, and, similar as in [n\ I18j. we reduce the arrangement computation to the geometric- 
topological analysis of a single curve and of a pair of curves. However, we improve in the following two 
main aspects: Firstly, we considerably reduce the amount of purely symbolic computations, that is, we 
only use resultant and gcd computation. The main reason for this approach is that we can outsource 
both computations to graphics hardware \20\ [21] |22] , removing a bottleneck of previous methods which 
was due to the high amount of symbolic operations. Secondly, for curve analysis, we use a result from 
Teissier |24[ I33j to obtain additional information for the number of distinct complex roots of f{a,y) along 
a critical fiber (actually, an upper bound which most likely matches the exact number). We combine 
this information with a new certified complex root solver |26| to isolate the roots of f{a, y). The latter 
symbolic-numeric step applies as an efficient filter denoted FastLift which fails only for very special 
instances. In order to achieve completeness of our overall method (i.e., for cases where FastLift fails), 
we modify the method from Uj for solving bivariate polynomial systems in order to isolate the roots of 
f{a,y)- 

We implemented our algorithm as a development branch of Cgal's^ bivariate algebraic kernel (Ak_2 for 
short) which is based on the algorithms from |17[ [T8] for topology and arrangement computation. Intensive 
benchmarks |18|. [29] have shown that Ak_2 can be considered as the current reference implementation. 
For fair comparison, we run an Ak_2 with GPU-enabled resultants and gcds against our implementation 
on numerous challenging benchmark instances. Our experiments show that the new approach outperforms 
Ak_2 for all instances. More precisely, our method is, on average, twice as fast for easy instances such 
as non-singular curves in generic position. For hard instances, we typically improve by large factors 
between 5 and 120 which is mainly due to the new symbolic- numeric filter FastLift, the exclusive usage of 
resultant and gcd as only symbolic operations and the abstinence of shearing. In summary, the presented 
approach demonstrates the strength of symbolic- numeric techniques. It further proves that, in order to 
achieve efficiency in an actual implementation, it is of great importance to reduce the amount of symbolic 
operations and to consider approximate operations whenever this is possible. 

Finally, we are confident that our new approach will have some positive impact in the following respects: 
Existing subdivision methods [2l lll [ [28 ] show excellent behavior for non-singular input. For singular input, 
they can be made certifying and complete when considering worst case separation bounds, however, this 
approach has not shown effective in practice so far. Another advantage of subdivision methods compared 
to elimination approaches is that they are local and do not need (global) algebraic operations. In our 
method, we considerably reduced the amount of such algebraic operations and outsourced the remaining 
computations to graphics card. Hence, it seems reasonable that combining our algorithm with a subdivision 
approach eventually leads to a certified and complete method which shows excellent "local" behavior as 
well. We further see numerous applications of our method, in particular, when computing arrangements of 
surfaces. The actual implementation ^ for surface triangulation is crucially based on planar arrangement 
computations of singular curves. Thus, we are confident to considerably improve its efficiency based on 
the new algorithm for planar arrangement computation. 



^Computational Geometry Algorithms Library, [www . cgal . org see also http://exacus.mpi-inf.mpg.de/cgi-bin/xalci 



|cgi for an online demo on arrangement computation. 
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Figure 2.1: The figure on the left shows a curve C with two x-extremal points and one singular point (red dots). In 
the projection phase, these points are projected onto the a;-axis and rational points separating the x-critical values 
are inserted (red crosses). In the lifting phase, the fibers at the critical values (red dots) and at the rational points in 
between (red crosses) are computed. In the connection phase, each pair of lifted points connected by an arc of C is 
determined and a corresponding line segment is inserted. The right figure shows the final graph that is isotopic to C. 

2 Curve Analysis 
2.1 The Algorithm 



The input of our algorithm is a planar algebraic curve C as defined in (1.1), where / is a square-free, 
bivariate polynomial / G Z[x,y] with integer coefficients. If / is considered as polynomial in y with 
coefficients fi{x) S its coefficients typically share a trivial content h := gcd(/o, /i, . . .). A non-trivial 
content h ^ Z defines vertical lines at the real roots of h. Our algorithm handles this situation by dividing 
out h first and finally merging the vertical lines defined by /i = and the curve analysis of the curve 
C := V{f /h). Hence, throughout the following considerations, we can assume that h is trivial, thus C 
contains no vertical line. 

The algorithm returns a planar graph Qc that is isotopic^ to C, where all vertices V of Qc are located 
on C. The proposed algorithm follows a classical cylindrical algebraic decomposition approach. We start 
with a high-level description of the three-step algorithm. 

In the first step, the projection phase, we project all x-critical points {a,P) G C (i.e., f{a,l3) = 
fy{a, (3) = 0) onto the x-axis by means of a resultant computation and root isolation for the elimination 
polynomial. The set of x-critical points comprises exactly the points where C has a vertical tangent or is 
singular. It is well known (e.g., see \25\ Theorem 2.2.10] for a short proof) that, for any two consecutive 
x-critical values a and a', C is delineable over / = [a, a'), that is, C|/xiR decomposes into a certain number 
mj of disjoint function graphs Cj^i, . . . ,Cj^rni- In the lifting phase, we first isolate the roots of the 
(square- free) intermediate polynomial f{qi,y) G Q[y], where qj constitutes an arbitrary chosen but fixed 
rational value in /. This computation yields the number nij (= number of real roots of f{qi,y)) of arcs 
above I and corresponding representatives {qi,yi,i) € Cj^i on each arc. We further compute all points 
on C that are located above an x-critical value a, that is, we determine the real roots ya,i, ■ ■ ■ ,ya,ma 
of each (non square-free) fiber polynomial f{a,y) E From the latter two computations, we obtain 

the vertex set V of Qc as the union of all points {qi,yi,i) and {a,ya^i). In the final connection phase, 
which concludes the topology analysis, we determine which of the above vertices are connected via an 
arc of C. For each connected pair (^1,^2) G V ■, we insert a line segment connecting vi and V2- It is then 
straightforward to prove that Qc is isotopic to C; see Theorem 6.4.4] for a proof. We remark that we 
never consider any kind of coordinate transformation, even in case where C contains two or more x-critical 
points sharing the same x-coordinate. 

We next describe the three phases in detail: 



^Qc is isotopic to C if there exists a continuous mapping : [0, 1] x C such that 0(0, C) — C, <j){l, C) = Qc and 

4>{to, .) : C i-> 4>{to, C) constitutes a homeomorphism for each to G [0, 1]. 



2.1.1 Projection Phase 

In the projection step, we follow well-known techniques from elimination theory, that is, we compute the 
resultant R{x) := ies{f,fy,y) £ 7L\x\ and a square-free factorization of R. More precisely, we determine 

square-free and pairwise coprime factors G i = 1, . . . ,deg(i?), such that R{x) = HS*'^^ (^i(^))*- 
We remark that, for some i G {1, . . . , deg(i?)}, rj(x) = 1. Yun's algorithm [35', Alg. 14.21] constructs 
such a square-free factorization by essentially computing greatest common divisors of R and its higher 
derivatives in an iterative way. Next, we isolate the real roots ajj, j = 1, . . . of the polynomials rj 
which in turn are i-fold roots of R. More precisely, we compute disjoint intervals I{aij) C M with rational 
endpoints such that I{ai,j) contains a^j but no other root of rj, and the union of all I{aij), j = 1, . . . 
contains all real roots of rj. For the real root isolation, we consider the Descartes method \15\ I30j as a 
suited algorithm. By further refining the isolating intervals, we can achieve that all intervals I{aij) are 
pairwise disjoint and, thus, also constitute isolating intervals for the real roots of R. Then, for each pair a 
and a' of consecutive roots of R defining an open interval / = (a, a'), we choose a separating rational 
value qj in between the corresponding isolating intervals. 

2.1.2 Lifting Phase 

Isolating the roots of the intermediate polynomials f{qi,y) is rather straightforward: since f{qi,y) is a 
square-free polynomial with rational coefhcients, the Descartes method directly applies. Determining the 
roots of f{a,y) S M[?/] at an x-critical value a is more complicated because f{a,y) has multiple roots and, 
in general, irrational coefficients. We propose to run the following approach: We first consider a method 
denoted FastLift which works as a filter for the fiber computation. We will see in the experiments enlisted 
in Section |4] that FastLift applies to all fibers for the majority of all input curves and only fails for a small 
number of fibers for some very special instances. In case of success, the fiber at a is returned. If FastLift 
fails, we use a second method denoted Lift which serves as a "backup" for FastLift. In comparison to 
FastLift, Lift is a complete method which applies to any input curve and any corresponding x-critical 
value, however, for the price of being less efficient. Nevertheless, our experiments show that even the 
exclusive use of Lift significantly improves upon existing approaches. 

Lift — a complete method for fiber computation: Lift is based on our recent studies on solving a 
bivariate polynomial system. In [4j, we introduced a highly efficient method, denoted Bisolve, to isolate 
the real solutions of a system of two bivariate polynomials f,g G Q[x,y]. Its output consists of a set of 
disjoint boxes Bi, . . . ,Bm C M? such that each box Bi contains exactly one real solution ^ := (xo,yo) 
of f{x,y) = g{x,y) = 0, and the union of all Bi covers all solutions. Furthermore, for each solution ^, 
Bisolve provides square-free polynomials p,q £ 7L\x\ with p(xo) = ^(yo) = and corresponding isolating 
(and refineable) intervals /(xq) and /(yo) for xq and respectively. Comparing ^ with another point 
'^1 = {xxiVi) £ given by a similar representation is rather straightforward. Namely, let p,q £ Z[x] 
be corresponding defining square-free polynomials and /(xi) and I{yi) isolating intervals for xi and yi, 
respectively, then we can compare the x- and y-coordinates of ^ and ^ via gcd-computation of the defining 
univariate polynomials and sign evaluation at the endpoints of the isolating intervals (see [3l Algorithm 
10.44] for more details). 

In order to compute the fiber at an x-critical value a of C, we proceed as follows: We first use Bisolve 
to determine all solutions pi = {a, /3j), i = 1, of the system f = fy = with x-coordinate a. Then, 
for each of these points, we compute 

ki := min{A; : fyk{a,^i) = ^("'/^O / 0} > 2. 

The latter computation is done by iteratively calling Bisolve for fy = fyi = 0, fy2 = fyS = 0, etc., and 
sorting the solutions along the vertical line x = a. We eventually obtain disjoint intervals /i, . . . , I; and 
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corresponding multiplicities ki,. . . ,ki such that f3j is a kj-iold root of f{a,y) which is contained in Ij. 
The intervals Ij already separate the roots f3j from any other multiple root of f{a,y), however, Ij might 
still contain ordinary roots of f{a,y). Hence, we further refine each Ij until we can guarantee via interval 

arithmetic that ^^-j^{a,y) does not vanish at Ij. If this condition is fulfilled, Ij cannot contain any root of 

f{a,y) except f3j due to the mean value theorem, thus, Ij is an isolating interval. 

After refining all intervals Ij, it remains to isolate the ordinary roots of f{a,y). For this purpose, we 
use the so-called Bitstream Descartes isolator |19j (Bdc for short) which can be considered as a variant 
of the Descartes method working on polynomials with interval coefficients. This method can be used to 
get arbitrary good approximations of the real roots of a polynomial with "bitstream" coefficients, that is, 
coefficients that can be approximated to arbitrary precision. Bdc starts from an interval guaranteed to 
contain all real roots of a polynomial and proceeds with interval subdivisions giving rise to a subdivision 
tree. Accordingly, the approximation precision for the coefficients is increased in each step of the algorithm. 
Each leaf of the tree is associated with an interval I and stores a lower bound I (I) and an upper bound 
u{I) on the number of real roots within this interval based on Descartes' Rule of Signs. Hence, u{I) = 
implies that / contains no root and thus can be discarded. If /(/) = u{I) = 1, then / is an isolating 
interval for a simple root. Intervals with u{I) > 1 are further subdivided. We remark that, after a number 
of iterations, Bdc isolates all simple roots of a bitstream polynomial, and intervals not containing any 
root are eventually discarded. For a multiple root ^, Bdc constructs intervals / which approximate ^ to 
an arbitrary good precision but never certifies such an interval / to be isolating. 

In our situation, we have already isolated the multiple roots of f{a, y) by the intervals Ij. It remains to 
isolate the simple roots of f{a,y). Therefor, we consider the following modification of Bdc : We discard 
an interval / if one of following three cases applies: i) u{I) = 0, or ii) I is completely contained in one 
of the intervals Ij, or iii) / contains an interval Ij and u{I) < kj. Namely, in each of these situations, / 
cannot contain an ordinary root of f{a,y). An interval / is stored as isolating for an ordinary root of 
/(a, y) if /(/) = u{I) = 1 and / intersects no interval Ij. All intervals which do not fulfill one of the above 
conditions are further subdivided. In a last step, we sort the intervals Ij (isolating the multiple roots) and 
the isolating intervals for the ordinary roots along the vertical line. 

We remark that Bisolve applied in Lift reuses the resultant obtained in the projection phase of the 
algorithm. Furthermore, it is a local approach in the sense that it suits very well for computing the fiber 
only at some specific j;-critical value a. Hence, when considering only a small number of x-critical fibers, 
the running time of Lift is considerably lower than its application to all fibers. 

FastLift — a fast method for fiber computation: FastLift is a hybrid method to isolate all complex roots 
and, thus, also the real roots of f{a,y), where a is an x-critical value of C. It combines a numerical solver 
to compute arbitrary good approximations (i.e., complex discs in C) of the roots, an exact certification step 
to certify the existence of roots within the computed discs and the following result due to Teissier \24:\ [33] : 

Lemma 1 (Teissier) For an x-critical point p = (a,/3) of C = V{f), it holds that 

mult(/(a,y),/3) = lnt{f,fy,p) - lnt{f^,fy,p) + 1, (2.1) 

where mult(/(a, y), /3) denotes the multiplicity of (3 as root of f{a,y) £ lnt{f,fy,p) the intersection 
multiplicity^ of the curves implicitly defined by f = and fy = Oat p, and Int(/a;, fy,p) the intersection 
multiplicity of fx = and fy = Oat p. 

Remark. In the case, where fx and fy share a common non-trivial factor h = gcd{ f x, fy) G 'L[x,y\, h 
does not vanish at any j;-critical point p of C. Namely, h{p) = would imply that lTit{fx, fy,p) = oo and, 

''The intersection multiplicity of two curves / = and p = at a point p is defined as the dimension of the localization of 
C[x,y]/ {f, g) at p, considered as a C-vector space. 
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thus, Int(/, fy,p) = oo as well, a contradiction to our assumption on / to be square-free. Hence, we have 
Int{fx, fy,p) = Int(/*,/*,p) with /* := fx/h and /* := fy/h and, thus, the following formula (which is 



equivalent to (2.1) for trivial h) applies: 

mult(/(a,y),/3) = Int(/,/^,p)-Int(/:,/;,p) + l, (2.2) 

We come to the description of FastLift. In the first step, we determine an upper bound m* for the 
actual number ma of distinct complex roots of f{a,y) which most likely matches m^. We distinguish the 
cases deg/(a, y) 7^ deg^^ / and degf{a,y) = degyf. In the first case, C has a vertical asymptote at a. 
Then, we define m* := deg/(a,y) which is obviously an upper bound for m^. If C is in generic position, 
f{a,y) has only ordinary roots and, thus, m* = m^. 



We now consider the case deg/(a, y) = deg„ /. Then, due to the formula (2.2), we have 



iTT-a = #{distinct complex roots of /(a, y)} 
= degyf - deg gcd{f {a, y)Jy{a,y)) 
= Yl (mult(/(a,y),/3)-l) 

/3eC:/(a,/3)=0 

= deg/,- (lnt(/,/„(a,/3))-Int(/:,/*,(a,/3))) 

/3GC: 
(a, (5) is x-critical 

= deg,/-mult(i?,a)+ ^ Int(/;, (a, (2.3) 

(a, 13) is x-critical 

< deg, / - mult(fl, a) + J2 Int(/:, /*, (a, /?)) (2.4) 

/3ec 

= degy f — mult(ii, a) + umlt{Q, a) =: m* (2-5) 



where R{x) = Tes{f,fy,y) and Q{x) := res(/*, /*, y). The equality (2.3) is due to the fact that / has no 
vertical asymptote at a and, thus, the multiplicity mult(iZ, a) equals the sum ^^^^ Int((/, /y, {a,/3)) of 
the intersection multiplicities of / and fy in the fiber at a. (2.5) follows by an analogous argument for 



the intersection multiplicities of /* and /* along the vertical line at a. From the square- free factorization 
of R, we already know mult(i?, a), and umlt{Q, a) can be determined by computing Q, its square-free 
factorization and checking whether a is a root of one of the factors. If the curve C is in generic position^. 



the inequality (|2.4|) becomes an equality because then /* and /* do not intersect in any point above a 



which is not located on C. Thus, in this case, we have ttt-q, = m* 



In the second step of FastLift, we aim to isolate all (complex) roots of /(a, y). The above computation 
of an upper bound m* motivates the following ansatz: In order to determine the roots of f{a,y), we 
combine a numerical complex solver for /(a, y) and an exact certification step. More precisely, we use the 
numerical solver to determine disjoint discs Di,. . . ,Dm in the complex space and an exact certification 
step to certify the existence of a certain number mj > 1 of roots (counted with multiplicity) of f{a, y) 



within each Df, see Section 2.2 for details. Increasing the working precision and the number of iterations 
for the numerical solver eventually leads to arbitrary well refined discs Di - but without a guarantee that 
these discs are actually isolating! Since m* constitutes an upper bound on the number of distinct complex 
roots of f(a,y), we must have m < < m* at any time. Hence, if the number of discs m equals the 
upper bound m* , we know for sure that all complex roots of deg/(a,y) are isolated. Then, the isolating 



^The reader may notice that generic position is used in a different context here. It is required that all intersection points 
of and fy above a are located on the curve C. We remark that this requirement is usually fulfilled even if there are two 
x-critical values sharing the same a::-coordinate. 
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a J = {a,a') a' a [{a) b I'ict) 



Figure 2.2: The left figure shows the generic case, where exactly one x-critical point (pa) above a exists. The 
bottom-up method connects Ai to pi and A2 to P2] the remaining arcs have to pass p^. In the second figure, the 
fiber at a contains two critical points p2 and p^. The red horizontal line segments pass through arbitrary chosen 
points (a,ti) separating pi-i and pi. The initial isolating interval I{a) = {a,b) for a is not sufficient to determine 
the connections for all arcs since ^1,^2,^3 intersect the segments / x {ti}. On the right, the refined isolating 
interval I' (a) induces boxes I' (a) x {ti,ti^i) small enough such that no arc crosses the horizontal boundaries. By 
examination of the y-coordinates of the intersections between the arcs and the fiber over the right-hand boundary of 
I' (a) (red crosses), we can match arcs and critical points. 



discs Di, . . . , Dm are further refined until, for alH = 1, . . . , m, 

A n M = or A n Dj = for all j / i, (2.6) 

where Di := {z : z G Di} denotes the complex conjugate of Di. The latter condition guarantees that each 
disc Di which intersects the real axis actually isolates a real root of f{a,y). In addition, for each real root 
isolated by some Di, we further obtain its multiplicity rrii as a root of /(a, y). 

If m ^ m* , either one of the roots of /(a, y) is still not isolated or it holds nia < rn*^. In the case that 
we do not succeed, that is, we still have m < m*^ after a number of iterations (in the implementation, this 
number is set empirically) with increasing precision, we stop FastLift and proceed with Lift instead. In 



Section 2.3 we discuss the few failure cases. 



We remark that the usage of Teissier's formula is not entirely new when computing the topology of 



algebraic curves. In |12| [29], the formula (2.1 ) was used in its simplified form to compute mult(/3, /(a, y)) 
for a non-singular point p = {a,j3) (i.e., lnt{fx, fy,p) = 0). In contrast, we use the formula in its general 
form and sum up the information along the entire fiber which eventually leads to the upper bound m* on 
the number of distinct complex roots of f{a,y). 



2.1.3 Connection Phase 

Let us consider a fixed x-critical value q, the corresponding isolating interval I{a) = (a, b) computed in 
the projection phase and the points pi := (q, ya,i) £ C, i = 1, . . . , nia, located on C above a. Furthermore, 
let / = {a, a') be the interval connecting a with the nearest x-critical value to the right of a (or -|-oo if 
none exists) and Aj, j = 1, . . . , mj, the j-th arc of C above / with respect to vertical ordering. To its left, 
Aj is either connected to (a, ±00) (in case of a vertical asymptote) or to one of the points pi. In order to 
determine the point to which an arc Aj is connected, we consider the following two distinct cases: 

• The generic case, in which there exists exactly one real x-critical point pj,, above a and deg f{a, y) = 
degy /. The latter condition implies that C has no vertical asymptote at a. Then, the points 
pi, . . . ,Pio^i must be connected with Ai, . . . ,^iQ_i in button-up fashion, respectively, since, for 
each of these points, there exists a single arc of C passing this point. The same argument shows 
that Pio+i, ■ ■ ■ ,Pma must be connected to Amj-ma+io+i^ • ■ ■ > in top-down fashion, respectively. 
Finally, the remaining arcs in between must all be connected to the x-critical point pjy. 

• The non-generic case: Each arc Aj is represented by a point Oj := {qi,yij) G C, where yjj denotes 
the j-th real root of f{qi,y) and qi an arbitrary but fixed rational value in I. We choose arbitrary 



7 



rational values h,..., tm^+i with ti < ya,i < t2 < . . . < Va^m^ < W+i- Then, the points pi := (a, U) 
separate the pj's from each other. Computing such pi is easy since we have isolating intervals with 
rational endpoints for each of the roots ya,i of f{a,y). In a second step, we use interval arithmetic 
to obtain intervals ^f{I{a) x ti) with f{I(a) x ti) C x ti). As long as there exists an i 

with E 5S/(/(a) x ti), we refine I (a). Since none of the pi is located on C, we eventually obtain a 
sufficiently refined interval I{a) with ^ x tj) for all i. It follows that none of the arcs 

Aj intersects any line segment I (a) x tj. Hence, above /(a), each Aj stays within the the rectangle 
bounded by the two segments 1(a) x ti^ and I{a) x tio+i and is thus connected to pi^. In order to 
determine zq, we compute the j-th real root 7^- of f{b, y) E Q[?/] and the largest zq such that 7^ > 
In the special case where 7^ < ti or 7^ > ti for all i, it follows that Aj is connected to {a, —00) or 
(a,+cxD), respectively. 

For the arcs located to the left of a, we proceed in exactly the same manner. This concludes the 
connection phase and, thus, the description of our algorithm. 

2.2 Numerical Solver with Certificate 

In FastLift, we deploy a certified numerical solver for a fiber polynomial to find regions certified to 
contain its complex roots. Bini and Fiorentino presented a highly efficient solution to this problem in their 
MPSOLVE package [8]. We adapt their approach in a way suited to handle the case where the coefficients are 
not known a priori, but rather in an intermediate representation which can be refined to any arbitrary finite 
precision. The description given in this section is high-level. For the details of an efficient implementation, 
we refer the reader to [26]. In the following considerations, let g{z) := f{a,z) = J27=o9i^^ ^ ^[-^1 be the 
fiber polynomial at an x-critical value a and V{g) = {Q}, i = 1, . . . , n, its complex roots. Thus, g{z) = 

nr=i(^-o). 

The main algorithm used in the numerical solver is the Aberth-Ehrlich iteration for simultaneous 
root finding. Starting from arbitrary distinct root guesses (zi)i=i^..,^„, it is given by the component-wise 
iteration rule z'^ = Zi if g{zi) = 0, and 

/ ^ 9{zi)/g'{zi) 

otherwise. As soon as the approximation vector lies in a sufficiently small neighborhood of some 
permutation of the actual roots {C,i)i of g, this iteration converges with cubic order ^34j- In practice, 
Aberth's method shows excellent performance even if started from an arbitrary configuration far away 
from the solutions. 

A straightforward implementation of [8| expects the coefficients gi of g to be known up to some 
relative precision p, that is, the input is a polynomial g = ^gix'^ whose floating point coefficients satisfy 
\gi — 9i| — l^il- In particular, this requirement implies that we have to decide in advance whether a 
coefficient vanishes. In general, though, the critical x-coordinate a of the fiber polynomial, and thus the 
coefficients of g, are not rational. Thus, it translates to expensive symbolic gcd computations of R and the 
coefficients of / as a univariate polynomial in Z[y][x]. 

Instead, we work on a Bitstream interval representation |16| [26] [gY of g. Its coefficients are interval 
approximations of the coefficients of g, where we require the width \gf —g^ \ of each coefficient [g]^ = [g^ ,g^] 
to be < /i for a certain absolute precision /i = 2~^. In this sense [g]^ represents the set {g '■ gi £ [g]i} of 
polynomials in a ^-polynomial neighborhood of g; in particular, g itself is contained in \g]^. Naturally, for 
the interval boundaries, we consider dyadic floating point numbers (bigfloats). Note that we can easily 
compute arbitrarily good Bitstream representations of /(a, z) by approximating a to an arbitrary small 
error, for example using the quadratic interval refinement technique [Ij. 

Starting with some precision (say, fi = 2~^^) and a vector of initial approximations, we perform 
Aberth's iteration on some representant g E [g]^- The natural choice is the median polynomial with 
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Figure 2.3: The two figures on the left show the topology analyses for the curves C — V{f) and D = V{g). The 
figure in the middle shows the intersection of the two curves. For the curve pair analysis, critical event lines (at dots) 
are sorted and non-critical event lines (at crosses) in between are inserted. Finally, for each event line a; = a, the 
roots of f{a,y) and g{a,y) are sorted. The latter task is done by further refining corresponding isolating intervals 
(blue or red intervals) and using the combinatorial information from the curve analyses and the computation of the 
intersection points. 

9i ~ + 9t)/'^^ take the liberty to select other candidates in case of numerical singularities in 

Aberth's rule (most notably, if g'{zi) = in some iteration). 

After a finite number of iterations (depending on the degree of (7), we interrupt the iteration and check 
whether the current approximation state already captures the structure of V{g). We use the following result 
by Neumaier and Rump [31], founded in the conceptually similar Weierstrafi-Durand-Kerner simultaneous 
root iteration: 



Lemma 2 (Neumaier) Let g{z) = gnY\^=i{z — Ci) G ^[z], g-n 7^ 0. Let Zi £ C fori = 
distinct root approximations. Then, all roots of g belong to the union D of the discs 



1. 



, n be pairwise 



Di := D{zi - ri, \ri\), 



where ri := — ■ — and : = 
2 gn 



9{zi) 



Moreover, every connected component CofD consisting of m discs contains exactly m zeros of g, counted 
with multiplicity. 

The above lemma applied to [g]'' using conservative interval arithmetic yields a superset C = {Ci, . . . , Cm} 
of regions and corresponding multiplicities Ai, . . . , Am such that, for each G C, all polynomials g G [g]^ 
(and, in particular, g) have exactly Afc roots in Ck counted with multiplicities. Furthermore, once the 
quality of the approximations {zi)i and [g]'^ is sufficiently high, C converges to V{g). 

In FastLift, where we aim to isolate the roots of g := f{a,y), we check whether m = m*. If the 
latter equality holds, Teissier's lemma guarantees that the regions Ck £ C are isolating for the roots of g, 
and we stop. Otherwise, we repeat Aberth's iteration after checking whether G [g]^^{zi). Informally, if 
this holds the quality of the root guess is not distinguishable from any (possibly better) guess within the 
current interval approximation of g, and we double the precision (/i' = /x^) for the next stage. 

Aberth's iteration lacks a proof for convergence in the general case and, thus, cannot be considered 
complete. However, we feel this is a merely theoretical issue: to the best of our knowledge, only artificially 
constructed, highly degenerate configurations of initial approximations render the algorithm to fail. 
Regardless of this assumption, the regions G C are certified to comprise the roots of g at any stage of 
the algorithm by Neumaier's lemma and the rigorous use of interval arithmetic. Furthermore, since we use 
the numerical method as a filter only, the completeness of the overall approach is not harmed. 



2.3 Discussion 

We conclude our description of the curve analysis algorithm with the following discussion on our method 
FastLift in the lifting step. FastLift is a certified method, that is, in case of success, it returns the 
mathematical correct result. However, the reader may notice that FastLift does not apply in all cases, a 
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reason why we additionally consider the complete but less efficient backup method Lift for some of the 
x-critical fibers. The failure of FastLift is either due to a very special geometric situation along a certain 
fiber or due to the behavior of the numerical solver. Special geometric situations are: 



(Geol) C has a vertical tangent at x = a and f{a,y) is not square-free, or 

fx „„J f* _ fv 

gCdU^Jy) '^"^ h - gCd(/,,/j,) 



(Geo2) there exists an intersection point of /* — a/i f \ and /* = 

„cdU f -1 above an x-critical value a 



of / which is not located on C . 

In case that none of the above special geometric situations is given (or removed by applying a shear; 
see below), the success of FastLift is guaranteed if the following conditions on the numerical solver are 
fulfilled: 

(Numl) The numerical solver is run for sufficiently many iterations, and 

(Num2) the approximations returned by the numerical solver converge against the roots of f{a,y) when 
increasing precision and number of iterations. 



As already mentioned in Section 2.2, we consider (Num2) as an exclusively theoretical problem. We 
further remark that, alternatively, we can use an exact and complete complex bitstream solver [32] 
to compute arbitrary good approximations of the fiber polynomial f{a,y), an approach which can be 
considered as an extension of Bdc to complex roots. However, for efficiency reasons, we decided to 
integrate a numerical method into our implementation instead. (Numl) is a practical problem and, in our 
implementation, we just empirically set the maximal number of iterations and the maximal precision as 
parameters depending on the degree and the bitsize of the polynomials. We noticed that the success of 
FastLift actually depends on these parameter values and consider it an interesting research question how 
they should be related to the given input to achieve optimal running times. 

It is worth mentioning that a complete and exact topology computation can be fully based on FastLift 
if we allow shearing (i.e., a coordinate transformation x i— )• x + for a s G Q). Namely, for all but except 
a finite number N of shearing factors,^ (Geol) and (Geo2) do not apply to the sheared curve. Hence, when 
considering + 1 pairwise distinct shearing factors one by one in circular order and increasing the number 
of iterations and the precision in the numerical solver for each factor, FastLift eventually succeeds for all 
fibers. 

In practice, as observed in our experiments presented in Section |4j the failure conditions for FastLift 
are almost negligible, as the method only fails for a few critical fibers on very special instances. For the 
remaining fibers and for all fibers of the majority of instances, FastLift is successful and extremely fast. 



3 Arrangement Computation 

Goal's recent implementation for computing arrangements of planar algebraic curves reduces all required 
geometric constructions (as intersections) and predicates (as comparisons of points and x-monotone curves) 
to the geometric-topological analysis of a single curve [18] and pairs of curves [T^; see also [5] and Goal's 
documentation [55] . 

Beyond the improved curve analysis proposed in Section [2] we aim to avoid subresultant sequences in 



general when analysing pairs of curves (see illustration in Figure 2.3), which is straightforward given the 
analyses of each single curve and the common intersection points of the two curves computed by Bisolve: 
Let C = V{f) and D = V{g) be two planar algebraic curves implicitly defined by the square-free 
polynomials /, g £ Z[x,y]. The curve analysis for C provides a set of critical event lines x = a, where 

^In [251 Corollary 3.2.10], it is shown that there exist at most n* + n shearing factors such that the sheared curve has 
covertical a;-critical points. In analogous manner, one can compute an upper bound for the number of shearing factors, where 
the sheared curve is not in (Geol) or (Geo2). 
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/(a, y) is non square-free. Each a is represented as the root of a square-free polynomial rj, with a factor 
of Rc '■= {f,fy,y), together with an isolating interval I (a). In addition, we have isolating intervals for 
the roots of f(a,y). A corresponding result also holds for the curve D with Rd := j:es{g , gy , y) . For the 
common intersection points of C and D, a similar representation is known. That is, we have critical event 
lines X = a', where a' is a root of a square-free factor of Rcd '■= ^es{f,g,y) and, thus, f{a,y) and g{a,y) 
share at least one common root (or the their leading coefficients both vanish for x = a). In addition, 
isolating intervals for each of these roots are computed. The curve-pair analysis now essentially follows 
from merging this information. More precisely, we first compute merged critical event lines (via sorting 
the roots of Rc, Rd and Rcd) and, then, pad merged non-critical event lines at rational values qj in 
between. The intersections of C and D with a non-critical event line at x = qj are easily computed by 
isolating the roots of f{qi,y) and g{qi,y) and further refining the isolating intervals until all isolating 
intervals are pairwise disjoint. For a critical event line x = a, we refine the already computed isolating 
intervals for f{a,y) and g{a,y) until the number of pairs of overlapping intervals matches the number m 
of intersection points of C and D above a. This number is obtained from the output of Bisolve applied 
to / and g, restricted to x = a. The information on how to connect the lifted points is provided by the 
curve analyses for C and D. 

We remark that, in the previous approach by Eigenwillig and Kerber [T7j, m is determined via efficient 
filter methods, too, while in general, a subresultant computation is needed if the filters fail. This is, for 
instance, the case when two covertical intersections of C and D occur. 



4 Implementation and Experiments 

Setup. We have implemented our algorithms as a branch of the bivariate algebraic kernel released with 
version 3.7 in October 2010, and replaced the curve and curve-pair analyses therein with our new methods 
based on FastLift, Lift and Bisolve. As throughout Cgal, we follow the generic programming paradigm, 
which allows us to choose among various number types and methods to isolated the real roots of integral 
univariate polynomials. For our setup, we rely on the number types provided by Gmp 5.0.1^ and the 
highly efficient univariate solver based on the Descartes method contained in Rs by Fabrice Rouillier |30j.^ 
which is also the basis for Isolate in Maple 13. 

All experiments have been conducted on 2.8 GHz 8-Core Intel Xeon W3530 with 8 MB of L2 cache on 
a Linux platform. For the GPU-part of the algorithm, we have used the GeForce GTX580 graphics card 
(Fermi Core). All used data sets are available for download.^ 

Symbolic Speedups. Our algorithm exclusively relies on two symbolic operations, that is, resultant 
and gcd computation. We outsource both computations to the graphics hardware to reduce the overhead 
of symbolic arithmetic which typically constitutes the main bottleneck in previous approaches. Besides a 
quick introduction given next, we refer the interested reader to [20^ [2T| l22j for more details. 

Shortly, both approaches are based on the "divide-conquer-combine" principle used in the modular 
algorithms by Brown [lOj and Collins [H] . This principle allows to distribute the computation over a large 
number of processor cores of the graphics card. At highest level, the modular approach can be formulated 
as follows. 1. Apply modular and/or evaluation homomorphisms to reduce the problem to a large set of 
subproblems over a simple domain. 2. Solve the subproblems individually in a finite field. 3. Recover the 
result with polynomial interpolation (in case evaluation homomorphism has been applied) and Chinese 
remaindering. 

Altogether, the GPU realization is quite straightforward and does not deserve much attention within 
the context of this work. It is only worth noting that our implementation of univariate gcds on the graphics 



Gm p: http : //gmplib . or g 



'Rs: 



http: //www. loria.f r/equipes/vegas /rs 



http: //www.mpi-inf .mpg. de/departments/di/projects/Geometry/DataSetsSNC-2011 .zip 
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card is comparable in speed with the one from Ntl^ 5.5 running on the host machine. Our intuition 
for this is because, in contrast to bivariate resultants, computing a gcd of moderate degree univariate 
polynomials does not provide a sufficient amount of parallelism, and Ntl's implementation is nearly 
optimal. Moreover, the time for initial modular reduction of polynomials, still performed on the CPU, 
can become noticeably large, thereby neglecting the efficiency of the GPU algorithm. Yet, we find it 
very promising to perform the modular reduction on the GPU which should further speed-up our algorithm. 

Contestants. We compare the new implementation with Cgal's bivariate algebraic kernel (see [5] 
and [6]) that has shown excellent performance in exhaustive experiments over existing approaches, namely 
CAD2d^^ and ISOTOP which is based on Rs. Both contestants were, except for few example instances, 
less efficient than Cgal's implementation, so that we omit further tests with them. Two further reasons 
can be given: Firstly, we enhanced Cgal's kernel with GPU-supported resultants and gcds, which makes 
it more competitive to existing software, but also to our new software, in case of non-singular curves, 
though slowdowns are still expected for singular curves or curves in non-generic position due to the need of 
subresultants sequences performed on the CPU. Secondly, the contestants based on Rs require as subtask 
Rs to solve the bivariate polynomial system / = /y = in the curve- analysis. In [3], we learned that the 
GPU-supported Bisolve was at least competitive to the current version of Rs and even showed in most 
cases an excellent speed gain over Rs. However, Rs is currently getting a very budding polish based on 
Rational Univariate Representations and modular arithmetic. Its theory will be presented at EuroCG 
2011 [9j. We are looking forward to compare our algorithms with a preliminary version of the new Rs 
quite soon. Moreover, we are confident that the authors will support the computation of arrangements of 
algebraic curves in a similar setup. Though tackling the same problem, both realizations can be seen as 
orthogonal if not complementary. 

4.1 Curve Analysis 

We first present the experiments comparing the analyses of single algebraic curves for different families of 
curves: (R) random curves of various degree and bit-lengths of their coefficients, (I) curves interpolated 
through points on a grid, (S) curves in the two-dimensional parameter space of a sphere, (T) curves that 
were constructed by multiplying a curve f{x,y) with f{x,y + 1), such that each fiber has more than 
one critical point. (P) projections of intersections of algebraic surfaces in 3D and, finally, (X) "special" 
curves of degrees up to 42 with many singularities or high-curvature points. We already considered these 
special curves in [1] and describe them in more detail in Table [s] in Appendix [A} The non-random and 
non-special curves are taken from [25', 4.3]. For the curve topology analysis, we consider four different 
setups: (a) Bisolve is, strictly speaking, not comparable with the actual curve analysis as it only computes 
the solutions of the system f = fy = 0. Still, it is interesting to see that our new hybrid method even 
outperforms this algorithm that only solves a subproblem of the curve- analysis, (b) LiftAna exclusively 
uses Lift for the fiber computations, (c) FastAna combines FastLift and Lift in the fiber computations. 
More precisely, it uses FastLift first, and if it fails for a certain fiber. Lift is considered for this fiber 
instead, (d) Ak_2 is the bivariate algebraic kernel shipped with Cgal 3.7, but with GPU-supported 
resultants and gcds. FastAna is our default setting, and its running time also includes the timing for the 
fiber computations where FastLift fails and Lift is applied instead. 

Table [T] lists the running times for single-curve analyses. We only give the results for representative 
examples; complete tables are given in Appendix [Aj It is easy to see that FastLift is generally superior to 
the existing kernel, even though Cgal's implementation profits from GPU-accelerated symbolic arithmetic. 
Moreover, while the speed-up for curves in generic position is already considerable (about half of the time), 
it becomes even more impressive for projected intersection curves of surfaces and "special" curves with 
singularities. The reason for this tremendous speed-up is that, for singular curves, Ak_2's performance 

I ^Ntl, http : //www . shoup . net/ntl^ 

http -.//mrw . usna . edu/Users/cs/qepcad/B/QEPCAD . html 
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Table 1: Running times (in sec) for analyses of algebraic curves of various families; timeout: algorithm timed out 
(> 400 sec) 



(R) sets of five random curves 


degree, bits 


BiSOLVE 


LiftAna 


FastAna 


Ak_2 


dense curves 


9, 10 


0.83 


0.73 


0.24 


0.66 


9, 2048 


4.17 


9.00 


2.24 


3.43 


15, 10 


2.82 


3.47 


1.01 


2.21 


15, 2048 


20.48 


50.06 


13.31 


16.82 


sparse curves 


9, 10 


0.25 


0.33 


0.11 


1.00 


9, 2048 


0.25 


0.39 


0.15 


0.98 


15, 10 


1.50 


3.31 


0.57 


3.19 


15, 2048 


15.66 


32.15 


8.08 


24.32 


(!) curve interpolated through points on a grid 


degree 


BiSOLVE 


LiftAna 


FastAna 


Ak_2 


9 


4.51 


7.25 


2.50 


4.98 


12 


25.25 


45.58 


14.66 


27.07 


(S) parametrized curve on a sphere 


with 16bit-coefficients 








degree 


BiSOLVE 


LiftAna 


FastAna 


Ak_2 


6 


6.52 


8.17 


2.49 


14.65 


9 


51.92 


78.09 


24.14 


45.66 


(T) curves wih a vertically translated copy 


degree 


BiSOLVE 


LiftAna 


FastAna 


Ak_2 


6 


5.27 


3.15 


0.83 


12.89 


9 


14.78 


14.91 


2.97 


159.22 


(P) projected intersection curve of surfaces with 8bit-coefficient 


degree(s) 


BiSOLVE 


LiftAna 


FastAna 


Ak_2 


2 ■ 2 


0.26 


0.23 


0.09 


0.22 


3 ■ 3 


2.87 


1.75 


0.52 


1.97 


4-4 


10.01 


11.93 


2.24 


38.91 


5 ■ 5 


83.00 


95.37 


13.76 


timeout 


(X) special curves (see Tablejs] AppendixjBjfor descriptions) 


name 


BiSOLVE 


LiftAna 


FastAna 


Ak_2 


curve_issac 


2.89 


2.30 


0.44 


2.96 


SA_2_4_eps 


0.35 


2.18 


0.64 


73.00 


grid_deg_10 


1.32 


2.52 


0.79 


1.54 


L6_circles 


4.52 


92.24 


2.10 


224.19 


huge_cusp 


8.50 


18.40 


5.50 


16.10 


swinnerton 


19.61 


16.81 


7.12 


304.09 


degree_7_surf 


14.57 


48.94 


7.39 


timeout 


spider 


57.09 


195.21 


26.26 


timeout 


FTT_5_4_4 


11.70 


44.27 


32.23 


timeout 


challenge_12b 


16.48 


52.71 


44.30 


timeout 


mignotte_xy 


260.61 


270.70 


136.54 


timeout 



drops significantly with the degree of the curve when the time to compute subresultants on the CPU 
becomes the dominating bottleneck of that approach. In addition, for curves in non-generic position, the 
efficiency of Ak_2 is affected because a coordinate transformation has to be considered in these cases. 

Recall that the symbolic-numerical filter in FastLift fails for very few instances, in which case 
Lift is locally used instead. The switch to the backup method is observable in timings; see for instance, 
challenge- 12b. As a result, the difference of the running times between Lift and FastLift are considerably 
less than for instances where the filter method succeeds for all fibers. In these cases, the numerical solver 



cannot isolate the roots within a given number of iterations, or we indeed have m < m* ; see Section 2.3 



Nevertheless, the running times are still very promising and yet perform better than Ak_2 for non- 
generic input, even though Lift's implementation is not yet mature enough, and we anticipate a further 
performance improvement. 
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Similar as Ak_2 has improved on previous approaches when it was presented in 2008, our new methods 
improve on Ak_2 now. That is, for random, interpolated and parametrized curves, the speed gain is 
noticeable, while for translated curves and projected intersections, we improve the more the higher the 
degrees. On special curves of large degree(!), we improve by a factors up to 100 and more. 



4.2 Arrangements 



For arrangements of algebraic curves, we compare two implementations: (A) Ak_2 is Cgal's bivariate 
algebraic kernel shipped with Cgal 3.7 but with GPU-supported resultants and gcds. (B) FastKernel 
is the same, but relies on FASxLiFT-filtered analyses of single algebraic curves. For the curve pair analysis, 
FastKernel exploits Ak_2's functionality whenever subresultant computations are not needed (i.e., a 
unique transversal intersection of two curves along a critical event line). For more difficult situations (i.e., 
two covertical intersections or a tangential intersection), the curve pair analysis uses BISOLVE as explained 
in Section [sj Our testbed consists of sets of curves from different families: (F) random rational functions of 
various degree (C) random circles (E) random ellipses (R) random curves of various degree and coefficient 
bit-length (P) sets of projected intersection curves of algebraic surfaces , and, finally, (X) combinations of 
"special" curves. 



We skip the tables for rational functions, 
circles, ellipses and random curves because the 
performance of both contestants are more or less 
equal: The linearly many curve- analyses are sim- 
ple and, for the quadratic number of curve-pair 
analyses, there are typically no multiple inter- 
sections along a fiber, that is, Bisolve is not 
triggered. Thus, the execution paths of both im- 
plementations are almost identical, but only as 
we enhanced Ak_2 with GPU-enabled resultants 
and gcds. In addition, we also do not expect 
the need of a shear for such curves, thus, the 
behavior is anticipated. The picture changes 
for projected intersection curves of surfaces and 
combinations of special curves whose running 
times are reported in Table [2} The Ak_2 re- 
quires for both sets expensive subresultants to 
analyze single curves and to compute covertical 
intersections, while FastKernel's performance 
is crucially less affected in such situations. 



Table 2: Running times (in sec) for computing arrangements of 
algebraic curves; timeout: algorithm timed out (> 4000 sec) 



(P) increasing number of projected surface intersections 



#resultants 


FastKernel 


Ak_2 


2 


0.21 


0.49 


3 


0.48 


0.93 


4 


1.03 


1.64 


5 


2.44 


3.92 


6 


5.14 


7.84 


7 


13.65 


21.70 


8 


22.69 


35.77 


9 


41.53 


67.00 


10 


58.37 


91.84 



(X) combinations of special curves (see Tablejs] Appendix B| 



#curves 


FastKernel 


Ak_2 


2 


9.2 


81.93 


3 


25.18 


148.45 


4 


248.87 


730.57 


5 


323.42 


836.43 


6 


689.39 


3030.27 


7 


757.94 


3313.27 


8 


1129.98 


timeout 


9 


1166.17 


timeout 


10 


1201.34 


timeout 


11 


2696.15 


timeout 
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A Further Experiments of Curve Analyses 



Table 3: Running times (in sec) for analyses of random algebraic curves; timeout: algorithm timed out (> 400 sec) 



(R) sets of five random curves 


degree, bits 


BiSOLVE 


LiftAna 


FastAna 


Ak_2 


dense curves 


6, 10 


0.38 


0.39 


0.15 


0.36 


6, 128 


0.32 


0.49 


0.18 


0.32 


6, 512 


0.55 


0.89 


0.35 


0.55 


6, 2048 


1.85 


3.50 


0.99 


1.71 


9, 10 


0.83 


0.73 


0.24 


0.65 


9, 128 


0.57 


0.90 


0.31 


0.55 


9, 512 


1.04 


1.91 


0.57 


0.96 


9, 2048 


4.17 


9.00 


2.24 


3.43 


12, 10 


2.18 


2.21 


0.69 


1.77 


12, 128 


1.64 


2.87 


0.86 


1.46 


12, 512 


2.84 


6.09 


1.56 


2.55 


12, 2048 


12.07 


30.26 


7.06 


9.96 


15, 10 


2.82 


3.47 


1.01 


2.21 


15, 128 


2.36 


4.55 


1.29 


2.06 


15, 512 


4.40 


9.84 


2.63 


3.77 


15, 2048 


20.48 


50.06 


13.31 


16.82 


sparse curves 


6, 10 


0.17 


0.21 


0.12 


0.22 


6, 128 


0.14 


0.19 


0.09 


0.22 


6, 512 


0.21 


0.36 


0.12 


0.39 


6, 2048 


0.73 


1.13 


0.39 


1.05 


9, 10 


0.25 


0.33 


0.11 


1.00 


9, 128 


0.25 


0.39 


0.15 


0.98 


9, 512 


0.43 


0.65 


0.23 


1.43 


9, 2048 


1.53 


2.51 


0.85 


4.31 


12, 10 


0.41 


0.83 


0.19 


1.67 


12, 128 


0.40 


1.05 


0.25 


1.60 


12, 512 


0.78 


2.16 


0.50 


2.53 


12, 2048 


3.82 


10.54 


2.82 


9.94 


15, 10 


1.50 


3.31 


0.57 


3.19 


15, 128 


1.48 


3.98 


0.72 


2.99 


15, 512 


2.96 


7.10 


1.45 


5.58 


15, 2048 


15.66 


32.15 


8.08 


24.32 
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Table 4: Running times (in sec) for analyses of algebraic curves of various families; timeout: algorithm timed out 
(> 400 sec) 



(1) curve interpolated through points on a grid 


degree 


BiSOLVE 


LiftAna 


FastAna 


Ak_2 


c 
O 


U.41 


U.oU 


U.zl 


n /in 
U.4y 


6 


0.72 


1.13 


0.41 


0.84 


■7 
1 


1 /lO 
l.^fZ 




n 7R 

U. f o 


1 fiO 
l.DZ 


3 


2.59 


3 78 


1.32 


2.86 


Q 


/I Ki 

't.Dl 


7 Ot^ 


Z.3U 


A OQ 


1 n 
lU 


D.Dz 


1 1 KO 


J.uz 


7 AC 
/ .Do 


11 




0"? 07 




1 j.yo 


Iz 




/IK E^Q 
4D.0O 


1 A fid 
14. DD 


07 n7 


13 


48.97 


93.97 


27.90 


46.21 


14 


101.96 


193.61 


59.14 


90.27 


15 


211.95 


timeout 


114.68 


166.68 


16 


timeout 


timeout 


236.39 


314.61 


(S) parametrized curve on a sphere 


with 16bit-coefFicients 








degree 


BiSOLVE 


LiftAna 


FastAna 


Ak_2 


6 


6.52 


8.17 


2.49 


14.65 


7 


22.58 


27.39 


8.47 


16.83 


8 


37.94 


53.3 


16.25 


32.09 


9 


51.92 


78.09 


24.14 


45.66 


10 


72.21 


110.81 


32.38 


64.31 


(T) curves with a vertically translated copy 


degree 


BiSOLVE 


LiftAna 


FastAna 


Ak_2 


5 


3.95 


1.93 


0.6 


5.61 


6 


5.27 


3.15 


0.83 


12.89 


7 


7.31 


6.14 


1.28 


33.97 


8 


9.66 


7.99 


1.75 


73.59 


9 


14.78 


14.91 


2.97 


159.22 


10 


17.03 


15.78 


3.36 


timeout 


(P) projected intersection curve of surfaces with Sbit-coefficient 


degree(s) 


BiSOLVE 


LiftAna 


FastAna 


Ak_2 


2 ■ 2 


0.26 


0.23 


0.09 


0.22 


2 ■ 3 


0.39 


0.40 


0.12 


0.32 


2 ■ 4 


1.03 


0.77 


0.25 


0.73 


2 ■ 5 


1.97 


1.53 


0.44 


1.93 


3 ■ 3 


2.87 


1.75 


0.52 


1.97 


3 ■ 4 


3.63 


3.02 


0.67 


5.21 


3 ■ 5 


8.71 


9.22 


1.90 


25.87 


4 ■ 4 


10.01 


11.93 


2.24 


38.91 


4-5 


21.45 


28.22 


4.46 


231.64 


5 ■ 5 


83.00 


95.37 


13.76 


timout 


(X) special curves (see Tablejs] AppendixjBjfor desciptions) 


name 


BiSOLVE 


LiftAna 


FastAna 


Ak_2 


curve.issac 


2.89 


2.30 


0.44 


2.96 


L4_circles 


1.56 


8.80 


0.52 


9.92 


SA_2_4_eps 


0.35 


2.18 


0.64 


73.00 


grid_deg_10 


1.32 


2.52 


0.79 


1.54 


ten.circles 


6.77 


3.79 


0.82 


24.07 


dfold_10_6 


4.45 


4.62 


1.58 


30.80 


L6_circles 


4.52 


92.24 


2.10 


224.19 


cov_sol_20 


8.68 


12.16 


3.38 


65.05 


curve24 


14.07 


15.61 


4.07 


50.73 


huge_cusp 


8.50 


18.40 


5.50 


16.10 


swinnerton 


19.61 


16.81 


7.12 


304.09 


degree_7_surf 


14.57 


48.94 


7.39 


timeout 


challenge_12a 


7.47 


15.84 


14.41 


timeout 


spider 


57.09 


195.21 


26.26 


timeout 


FTT_5_4_4 


11.70 


44.27 


32.23 


timeout 


challenge_12b 


16.48 


52.71 


44.30 


timeout 


mignotte_xy 


260.61 


270.70 


136.54 


timeout 
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B Description of Special Curves 



Table 5: Description of the special curves used in the experiments of Section |4j Sources of curves are given where 
known. 



Ills La 1 ]<_c 




L/ csi^l 1 p LltJI 1 


*_UI VC_I33C1*_ 


15 


ISLIIclLCLI pLIIIILS, 1 1 Igl 1 1 vd LU 1 C UUIIIL3 | ± 


L4_circles 


1 ft 

ID 


4 circles w.r.t. L4-norm| 






clustered solutions 


SA_2_4_eps 


16 


singular points with high tangencies, displaced [27| 


grid_dGg_10 




large coefficients! 






curve in generic position 


ten .circles 




set of 10 random circles multiplied together; rational solutions 


aroiu-iU-D 


^n 


many half-branches [27 


L6 -Circles 




4 circles w.r.t. L6-norm| 






cl ustered solutions 




20 


l—UVCl LIL.id 1 ^KJl U LltJI Is 


curve24 


24 


curvature of degree 8 curve; 






many singularities 


huge.cusp 


8 


large coefficients; 






high-curvature points 


swinnerton 


25 


covertical solutions in x and y 


degree_7_surf 


42 


silhouette of an algebraic surface; covertical solutions in x and y 


challenge_12a 


30 


many candidate solutions to be checked |27( 


spider 


12 


degenerate curve; 






many clustered solutions 


FTT_5_4_4 


40 


many non-rational singularities |27| 


challenge_12b 


40 


many candidates to check |27| 


mignotte_xy 


42 


a product of x'/i^-Mignotte polynomials, displaced; 






many clustered solutions 
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